Exact low-temperature behavior of kagome antiferromagnet at high fields 
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Low-energy degrees of freedom of a spin- 1/2 kagome antiferromagnet in the vicinity of the satura- 
tion field are mapped to a hard-hexagon model on a triangular lattice. The latter model is exactly 
solvable. The presented mapping allows to obtain quantitative description of the magnetothermo- 
dynamics of a quantum kagome antiferromagnet up to exponentially small corrections as well as 
predict the critical behavior for the transition into a magnon crystal state. Analogous mapping is 
presented for the sawtooth chain, which is mapped onto a model of classical hard dimers on a chain. 
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In the past few years there have been increased at- 
tention to the behavior of frustrated magnetic materials 
at high fields. The interest is stimulated by a search 
for new fundamental effects such as the magnetization 
plateauSfi*2i^*^iS*SiLSiS the magnetization jumps/'^-^-'^'-'^^ii^ 
new complicated magnetic phase diagrams and new 
types of the critical behavior^^ as well as by a pos- 
sible technological application of frustrated magnets 
in the adiabatic demagnetization refrigerators due to 
their unique magnetocaloric propertiesiS*ii A zero- 
temperature behavior of geometrically frustrated anti- 
ferromagnets in the vicinity of the saturation field He 
has been recently discussed by Schulenburg and co- 
workersiiSiiiiiSii^ Above He in the fully polarized phase, 
these models have a flat branch of magnons with energy 
{H — He). Such dispersionless excitations correspond in 
real space to localized spin flips, which do not interact 
with each other unless touched. Upon decreasing mag- 
netic field through He, the energy of a single excitation 
becomes negative and magnons condense into a close- 
packed crystal structure. The present work is devoted to 
theoretical investigation of the low-temperature thermo- 
dynamics of magnon crystals. Specifically, we consider 
spin- 1/2 Heisenberg models on a kagome lattice and a 
sawtooth (or A) chain, see Fig. ^ While the former lat- 
tice is one of the best known examples of geometric frus- 
tration, which is realized, e.g., in SrCr9pGai2-9pOi9 
the latter spin model is applicable to magnetic delafos- 
site YCu02.5!^ We show that the thermodynamic po- 
tential in the vicinity of He for the two quantum models 
can be calculated exactly up to corrections, which are 
exponentially small at low temperatures. 

We consider nearest- neighbor spin- 1/2 Heisenberg 
models in an external field 



excitations in the fully polarized phase is 
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For a kagome lattice all nearest-neighbor bonds have the 
same strength Jij = 1. The spectrum of single-magnon 



t^ik 



H-3, c^2,3k = - 1 ± 5^3(1 + 27k) , (2) 



where 7k = ^ 



is a sum over nearest-neighbor 



sites on a triangular Bravais lattice. At the saturation 
field He = i the energy of excitations from the lowest 
dispersionless branch vanishes. The sawtooth chain ex- 
hibits a similar behavior for the special ratio of the two 
coupling constants J2 — \Ji- In the following we always 
assume the above choice of the parameters with Ji = 1. 
Above He — 2 the sawtooth chain antiferromagnet has 
two branches of single-magnon excitations: 
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Excitations in the dispersionless branches correspond to 
localized states. Simplest examples of such states are 
spin flips trapped on hexagon voids of a kagome lattice 
or in the valleys of the sawtooth chain, see Fig. ^ Their 
wave functions are given by 



= 7sE(-l)""'^nJFM), 



(4) 



for the kagome and the sawtooth models, respectively, 
|FM) being the ferromagnetic state. Other more ex- 
tended localized states can be represented as linear com- 
binations of the above elementary states. Moreover, for a 
kagome lattice with N sites, there is the same number of 
hexagons as the number of states in the flat branch. 
Therefore, the localized states on the hexagons states 
can be used as a real-space basis for the lowest branch 
in Eq. ||2J). Similar correspondence exists between the 
\N states in the dispersionless branch Eq. JSjl and the 
localized states in the valleys of the sawtooth chain. 

Localized one-magnon states allow to construct a class 
of exact eigenstates in every n-magnon subsector with 
n < A^max by putting localized spin flips on isolated 
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FIG. 1: (color online) Kagome lattice (top) and the saw- 
tooth chain (bottom) with localized magnons. Amplitudes 
and phases of spin-down states are indicated near each site. 

hexagons or valleys. The highest possible density of inde- 
pendent localized magnons on the kagome lattice is one- 
third of the total number of hexagons or TVmax = ■^N. 
Such a magnon crystal state is three-fold degenerate and 
breaks the translational symmetry. A magnon crystal 
for the sawtooth chain has iVmax — and is two- fold 
degenerate. Since isolated localized magnons do not in- 
teract with each other, it is natural that they correspond 
to the lowest energy states in every n-magnon subsector 
(n < A^max)' "^'^'^^'^^'^^ Our idea is that only isolated lo- 
calized magnons contribute to the low temperature ther- 
modynamics. In order to prove this we have to show 
that (i) there are no other low energy states and that 
(ii) isolated localized magnon states are separated by a 
finite gap from the higher-energy states in the same S^- 
subsector. In the following we first prove the above two 
points and then show that localized magnons are mapped 
to statistical mechanics models of hard classical objects, 
hexagons or dimers, which allow exact solution for the 
thermodynamic properties. 

Numerical evidences and analytical arguments that 
isolated localized magnon states are the lowest energy 
states in every n-magnon sector have been presented 
beforeiifliiiiiSii^ Below we give a shortened proof for the 
spin- 1/2 kagome antiferromagnet, which is also valid for 
the checkerboard and the pyrochlore lattice antiferro- 
magnets. The Heisenberg Hamiltonian ^ is split in a 
standard way into the Ising part 7i^^, which includes 
all S^S^ terms plus the Zeeman energy, and the trans- 
verse part H'^. The energy of the ferromagnetic state 
is always subtracted from Ti.^^. A quantum state of n 
isolated localized magnons is not only an eigenstate of H 
with En = Ti{H — 3), but is also an eigenstate of 7i^^ and 
n-^ separately with E^^ = n{H - 2) and E^ = -n. The 



idea is to show that for an arbitrary state (7i^^) > E^^ 
and {H'^) > E:^. For Ti.^^ the formulated relation holds 
trivially: once all spin fiips sit on different bonds the Ising 
part of the Hamiltonian has the minimal energy E^^ . In 
order to show that the same is true for the transverse 
part we map the subspace of n spin flips on the Hilbert 
space of n hard-core bosons S„. The transverse part of 
the Heisenberg Hamiltonian is, then, the kinetic energy of 
bosons Ti.-^ = K,. In addition, we define the Hilbert space 
of n bosons without the hard-core constraint: B{)n 13 Bn- 
The minimum of the kinetic energy is easily found in Bon'- 
all n particles should occupy one of the low-energy single- 
particle states given by Eq. If the allowed quantum 
states are restricted again to i3„, from the variational 
principle the expectation value of K. can only increase: 
(^^) > ™inBon(^) = ^n- The presented arguments not 
only show that isolated localized magnons have the low- 
est energy, but they also allow to check whether these 
states are the only lowest energy states (see below). 

For the sawtooth chain the two branches of excita- 
tions in Eq. (jSJ are separated by a gap. The low-energy 
scattering states are, therefore, formed by magnons 
from the fiat branch. In the two-magnon sector there 
are ^{N/2){N/2 + 1) states constructed from N/2 one- 
magnon states of the lowest branch. The number of iso- 
lated localized magnons in this sector is ^{N/2){N/2 — 3). 
The difference gives the number of scattering states, 
which is equal to A^. In real space representation the 
above scattering states correspond to localized magnons 
occupying the same valley or adjacent valleys. This rep- 
resentation suggests the following basis in the subspace of 
scattering states: {ipi) = \(piipn^i) and jV'i) = l^pf), where 
the exclusion principle of spin-flips and proper normal- 
ization are implied. These states have nonzero overlaps: 

(V'»lV'«±i> = 1/35, = (V'»|V^»+i> - V5749. (5) 

We calculate the matrix elements of the Hamiltonian be- 
tween these states and use the variational principle to 
determine the lowest possible energy. A simple varia- 
tional state |fc) = lias lli^ energy E{k) = 
(fc|7i|fc)/(fc|fc). The minimal value of E{k) is reached at 
fc = with the gap A = |y w 0.54, which separates 
the two-magnon scattering states from the low-energy 
boundary 2{H — 2). An improved variational ansatz 

|fc) = Eie"*'"''(|V'») + c|'0»)) yields A « 0.44, which com- 
pares well with the exact diagonalization result 0.42i^ 
The gap is determined mostly by a nearest-neighbor re- 
pulsion of localized magnons in adjacent valleys. The 
dispersion of bound two-magnon complexes E(k) is weak 
and does not exceed 10%. We conjecture that the re- 
pulsion leads to a similar behavior in all n-magnon sec- 
tors: the lowest energy states with _E„ = n{H — 2) 
are separated from the scattering states by a finite gap 
A„ = 0(1). At temperatures T <C A„ one can neglect 
the higher energy states and consider only contribution of 
isolated localized magnons. The latter problem is equiv- 
alent to a one-dimensional lattice gas of particles with 
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energies {H — 2) and on-site and nearest-neighbor exclu- 
sion principle (classical hard-dimer model). 

The above consideration is straightforwardly extended 
to the kagome antiferromagnet. The two-magnon sec- 
tor consists of ^{N/3){N/3 + 1) states constructed from 
one-magnon states of the flat branch o^ik- They include 
^{N/3){N/3 — 7) states constructed by placing local- 
ized spin flips on isolated hexagons (called below hard- 
hexagon states) and the additional | A'' defect states. The 
latter states correspond to N states of localized magnons, 
which occupy adjacent hexagons, and to states with 
two spin flips on the same hexagon. We define the real- 
space basis for the defect states as \ipai) = \ViVi+Sc) 
and \^^) = \^^), where 6i = (1,0) and ,52,3 - (i±^) 
are three nearest-neighbor sites on a triangular lattice 
of hexagons. Using a simple variational ansatz |k) = 



-ikri 



Cal'tpai) we find the gap A « 0.24 between the 



isolated localized states and the scattering defect states. 

The low-energy sector of the kagome antiferromagnet 
has also two additional features. First, the hard-hexagon 
states are not the only independent localized magnon 
states. An extra two-magnon state is illustrated in the 
left part of Fig. ^ and consists of two spin flips cycling 
around a small and a large hexagons. Such a state is a 
special combination of the defect states \ijjai) and 
Combinatorial arguments show that the extra states give 
only corrections to the hard- hexagon result in n- 
magnon sectors with small n and disappear completely 
in sectors with 3n/N > 0.1 

Second difference with the sawtooth chain is that the 
gap for one of the dispersive branches in Eq. ^ also 
vanishes at H = He- The corresponding propagating 
magnon has the same energy wsk^o = -ff — 3 as local- 
ized magnons from the flat branch. The ground states 
in the two-magnon sector contain apart from the iso- 
lated localized magnons also superposition states of one 
localized magnon and one propagating magnon. Such 
states form a continuum above the low-energy thresh- 
old E2 = 2{H — 3). The same is true for all n-magnon 
sectors with n <^ N . Once n becomes a finite fraction 
of A'^, i.e. for partially magnetized states, the above pic- 
ture changes. The low-energy propagating magnons ex- 
perience multiple scattering from an infinite number of 
localized magnons. Such an interaction produces a fi- 
nite shift of their energy. The exact value of the en- 
ergy shift for the depends on a precise pattern of 
localized magnons. It can be explicitly estimated for 
a 'uniform' state of n-localized magnons, correspond- 
ing to a low-temperature ensemble of states with differ- 
ent translational patterns. We bosonize the spin Hamil- 
tonian using, e.g., the Holstein-Primakoff transforma- 
tions. The quadratic terms give the excitation spec- 
trum, which coincides with Eq. ||2Jl. To treat the ef- 
fect of interactions we define two Hartree-Fock averages 
TOi = {b\bi) and TO2 = (blbj). After decoupling the four- 
boson terms the effect of interaction is reduced to the 
renormalization S" — s- (S* — mi + TO2) in the quadratic 
terms, where S = 1/2 is a local spin. At the saturation 



field He = GS = 3 the gap for the dispersive mode is 
Ad = 6 (mi — 7712). For a single localized magnon one 
finds nil = — m2 = ^. Therefore, the dispersive mode 
opens a finite gap — 3n/N. At low temperatures we 
neglect these higher-energy states and reduce the isolated 
localized magnons on small hexagons to a gas of particles 
on a triangular lattice with on-site and nearest-neighbor 
exclusion, which is also called a hard-hexagon model. A 
remarkable feature of this model is that it has an ex- 
act solution, which was obtained by BaxterfSi Below, we 
discuss the low temperature properties of the kagome an- 
tiferromagnet and the sawtooth chain in the framework 
of the two exactly solvable models. 

The partition function of the classical hard dimers is 



^exp 



T 



(6) 



where e ~ H — He is the chemical potential and (Ti = 0, 1 
are the occupation numbers for dimers the index i runs 
over the N/2 sites of the basal chain. The nearest- 
neighbor exclusion principle is imposed by the last term. 
The partition function lO can be rewritten in terms of 
the transfer matrices T{<Ji, Ci+i ) as Z = TrT^/2. In the 
thermodynamic limit N 00 the free energy is deter- 
mined by the largest eigenvalue of the transfer matrix: 



F/N 
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(7) 



From this expression one can calculate various thermo- 
dynamic properties of the sawtooth chain. At H = He 
the entropy and the magnetization have universal values 
independent of temperature: 



M/N.i^. ,8, 



The entropy as a function of field reaches a sharp max- 
imum at H = He, which amounts to 34.7% of the to- 
tal entropy A^ln2 of the sawtooth chain. At low tem- 
peratures all the magnetization curves M{H) cross at 
H = He- Comparison between analytical and numeri- 
cal exact diagonalization results shows a nice agreement 
below T* « 0.1,^° The correlation length of the hard 
dimers is obtained from the ratio of the two eigenvalues 
of T, and the result is 



Vl + 46"-'/'^ + 1 



VT+4F^- r 



(9) 



with the wavenumber tt. 

The partition function of the hard-hexagon model is 
given by the same expression (|SJ), where the site in- 
dex runs over triangular lattice formed by the centers 
of hexagon voids of the original kagome lattice. This 
model was solved by Baxter, who used a corner transfer 
matrix method to obtain Z as a function of the fugacity 
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FIG. 2: Field dependence of the specific heat per one site (in 
units of fes) of a spin-1/2 kagome antiferromagnet at T = 0.05 
obtained in the hard-hexagon model. The inset shows the 
magnetization at the same T. 

z = e^^l"^ ?^ A number of exact results following from 
the Baxter's solution can be transferred to the spin-1/2 
kagome antiferromagnet using a trivial rescaling. The 
entropy as a function of field at constant T reaches a 
maximum at H = He- At the saturation field both the 
entropy and the magnetization have temperature inde- 
pendent universal values: 

S/N = 0.11108 , M/N ^ 0.44596 . (10) 

More important results concern a phase transition, which 
corresponds to a weak crystallization of hard hexagons or 
localized magnons in the two models: 

Hc{T) = Ha-Tlnzc , = (H + 5\/5)/2 , 
M,/N = i - ip, , = (5 - V5)/10 . (11) 

The second order transition at Hc{T) belongs to the uni- 
versality class of the three-state Potts model and has the 



exact critical exponents a = |, /? = |, and ?/ = |. 
The transition field has a linear temperature depen- 
dence, while the magnetization Mc stays along the crit- 
ical line. Variations of the specific heat and the magne- 
tization of the spin-1/2 kagome antiferromagnet derived 
from the hard-hexagon model^ are presented in Fig. El 
Baxter used the order parameter defined as a difference 
of hexagon densities pi on adjacent sites: R = pi — pi+s- 
Up to a renormalization prefactor R is related to the 
Fourier harmonics pq at q = (47r/3a*,0), where o* is a 
period of the triangular lattice. The magnon crystal has 
two types of spins, on hexagons with localized magnons 
and between them, which have different average mag- 
netizations: (S^) = 1/3 and 1/2 at T = 0, respectively. 
Hence, the crystalline order of magnons is reflected in ap- 
pearance of the corresponding Fourier harmonics in the 
structure factor S'^^(q) and can be in principle observed 
by elastic neutron scattering. 

In conclusion, the above analysis of the magnetother- 
modynamics of the quantum kagome antiferromagnet 
and of the sawtooth chain model has shown that strongly 
correlated spin systems have emergent "paramagnetic" 
degrees of freedom, which experience an effective mag- 
netic field h = H — He- Such effective "paramag- 
netic moments" are responsible for an enhanced magne- 
tocaloric effect of quantum geometrically frustrated an- 
tiferromagnets in the vicinity of the saturation field very 
similar to their classical counterparts^ During prepa- 
ration of the manuscript, we have learned about a re- 
cent preprint, where the authors have obtained a zero- 
temperature entropy of the two models in agreement with 
our Eqs. ^ and H10|l . though they have not discussed the 
low-temperature thermodynamics. 
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